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Abstract 

When the gauge fields have derivative couplings to scalars, like in the case of the rel¬ 
ativistic theory of Van der Waals (or Casimir-Polder) interactions, conformal invariance is 
broken but the magnetic and electric susceptibilities are not bound to coincide. We ana¬ 
lyze the formation of large-scale magnetic fields in slow-roll inflation and find that they are 
generated at the level of a few hundredths of a nG and over typical length scales between 
few Mpc and 100 Mpc. Using a new time parametrization that reduces to conformal time 
but only for coincident susceptibilities, the gauge action is quantized while the evolution 
equations of the corresponding mode functions are more easily solvable. The power spectra 
depend on the normalized rates of variation of the two susceptibilities (or of the correspond¬ 
ing gauge couplings) and on the absolute value of their ratio at the beginning of inflation. 
We pin down explicit regions in the parameter space where all the physical requirements (i.e. 
the backreaction constraints, the magnetogenesis bounds and the naturalness of the initial 
conditions of the scenario) are jointly satisfied. Weakly coupled initial data are favoured if 
the gauge couplings are of the same order at the end of inflation. Duality is systematically 
used to simplify the analysis of the wide parameter space of the model. 
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1 Introduction 


Magnetic fields with typical correlation scales exceeding the astronomical unit (i.e. roughly 
10 13 cm) permeate the interstellar and intergalactic plasmas which are, in many respects, very 
similar to the one we can produce in terrestrial experiments DEI- In spite of this, the origin 
of large-scale magnetism is still under intense debate both theoretically and observationally 
[3]. As proposed in the last few years [3] the temperature and the polarization anisotropies 
of the Cosmic Microwave Background (CMB in what follows) may be magnetized. This 
observation offers the unique opportunity of direct limits on the large-scale magnetism prior 
to matter-radiation equality since the large-scale magnetic fields affect directly the initial 
conditions of the Einstein-Boltzmann hierarchy. The current Planck explorer data can be 
used to set bounds on large-scale magnetic fields |5] as previoulsy done with the WMAP 
3-yr and 9-yr releases (see, respectively, first and second papers of Ref. m- The results can 
be summarized by saying that the WMAP9 Ed and Planck data [8] are compatible, they 
are both sensitive to magnetic fields in the nG range (1 nG = ICR 9 G) for magnetic spectral 
indices rig = 0(1.3) using the available temperature and polarization power spectra^ 

It has been repeatedly argued that large-scale magnetic fields might well be generated 
during a stage of inflationary expansion. The rationale for this requirement has to do with 
the correlation scale of the produced field that must be sufficiently large and the onset of the 
rotation of the protogalaxy. A very promising framework for generating magnetic fields with 
comoving correlation scales exceeding the Mpc is represented by models where the gauge 
fields couple directly to one or more scalar fields (see j9j M, HI H2] for an incomplete list of 
references). The scalar fields may coincide with one or (more inflatons) or even with multiple 
spectator fields. The conventional class of models is based on the following action)^} 


5 = 


167T 


d A x y/—g 


\{pA) Y Q pY a0 + ^■pseudoppt 'P')^a/3 V ^ 


( 1 . 1 ) 


where (p and ip may denote, for instance, a generic inflaton field and a generic spectator 
field. Various situations can be envisaged and most of them have been investigated in the 
literature. The presence of A (pj'ip) in Eq. (1.1) is more relevant than the pseudoscalar 
(axion-like US E]) coupling which will be ignored even if it has been studied by many 
authors pa ESI HZ! in the context of the magnetic field generation. For the amplification 
of the magnetic field itself the pseudoscalar vertex is not so efficient but it is relevant when 
one wants to generate magnetic fields whose flux lines are linked or twisted as originally 


2 We use here the same conventions employed to assign the curvature power spectra: the scale invariant 
magnetic power spectra are realized for tib —> i mm. Within these conventions the magnetic power 
spectrum (i.e. the Fourier transform of the two-point function) has the same dimensions of the magnetic 
energy density. While other conventions stipulate that the scale-invariant limit is realized for ns —> —3, in 
the present paper, the Fourier transform will be consistently assigned for the scalar modes of the geometry 
and for the magnetic fields implying that ub —> 1 is the scale-invariant limit. 

3 In the notations employed on this paper Y I ‘ L ' and Y^ v are, respectively, the gauge field strength and its 
dual; g = det g^ is the determinant of the four-dimensional metric with signature mostly minus. 
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discussed in HBj. The produced Chern-Simons condensate leads to a viable mechanism 
for baryogenesis via hypermagnetic knots [18j 19]. These helical fields play also a role 
in anomalous magnetohydrodynamics where the evolution of the magnetic fields at finite 
conductivity is analyzed in the presence of anomalous charges [20] , In the collisions of heavy 
ions this phenomenon is often dubbed chiral magnetic effect [2T |. 

It has been recently argued [22] that the class of models pinned down by Eq. (1.1) can 
be complemented by further terms: 


S = 


167r 


J d 4 x yf—g 


M^^)Y pa Y° 


NZM)Y pa Y° a 


( 1 . 2 ) 


Equation (1.2) leads to unequal electric and magnetic susceptibilities [22]. As a special 
case Eq. (1.2) includes the typical derivative coupling arising in the relativistic theory of 


Casimir-Polder and Van der Waals interactions 


S = 


I d'x V=g Ld aV >d eV >‘ Y°" Y^ + g 2 \ V \ 2 Y a „ Y^ 


(1.3) 


Other terms potentially present in Eq. (1.2) (such as V p Y aa Y pa ) will be neglected even if, 
as explained above, they might be relevant for the evolution of the magnetic helicity. 

The simplest parametrization of the coupling functions M. pa and M pa is [22] : 


M pa (ip) = \ E (ip)u p (ip)u a ((p), Afpa^) = Xsi^Up^Ua^). (1.4) 


Note that u p ((p) and u p {tJj) appearing in Eq. (1.4) are the normalized gradients of the 
corresponding scalar fields but, in the context of a purely hydrodynamical model, they can 
also play the role of the four-velocities of a relativistic fluid. The full action has been taken 
to be the sum of the first term of Eq. (1.1) and of the remaining two terms in Eq. (1.3) with 
the parametrization of Eq. (1.4). In this case the electric and the magnetic susceptibilities 
are given by: 

/ \ / \ 

(1.5) 




In this framework viable magnetogeneis models can be formulated in different dynamical 
situations and some of these possibilities have been already swiftly examined [22]. The 
purpose of the present paper is to undertake a comprehensive analysis of the parameter 
space of this scenario in terms of the initial conditions and of the evolution of the gauge 
couplings. To achieve this goal the expressions of A4 pa and M pa shall be first generalized. 
The three essential parameters of this class of models are the normalized rate of variation of 
the electric and magnetic gauge couplings (denoted, respectively, by Fe and Fb) and also 
the initial value of their ratio. 

From a pragmatic viewpoint, the whole parameter space can be reduced to the first 
quadrant of the ( F B , F E ) plane where F B and F E are both positive: the other three quadrants 
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can be charted through the systematic use of the duality transformation^] To save time we 
first impose the backreaction constraints over different scales, then discuss the naturalness of 
the initial conditions and finally verify if the obtained region of the parameter space satisfies 
the magnetogensis requirements. 

A technical aspect of the present analysis concerns the quantization of the action and the 
evolution of the related mode functions. This problem can be discussed in the conventional 
manner by using the conformal time coordinate r. However the analysis becomes much more 
transparent in terms of a newly time variable, denoted hereunder by r/, which reduces to the 
conformal time but only in the case of coincident electric and the magnetic gauge couplings. 

This paper is organized as follows. In section [2] we present the generally covariant decom¬ 
position of the coupling functions Ai pa and J\f prT and study their symmetry properties in the 
light of the equations of motion. We also show how the obtained parametrization fits within 
large held and small held inhationary models. In section [3] we show how the power spectra 
can be derived by avoiding the standard conformal time parametrization and by expressing 
the full action in terms of a new time variable directly related to the mismatch between 
the magnetic and the electric susceptibilities. The magnetic and electric power spectra are 
derived in section |4j we shall also analyze how the different regions of the parameter space 
are transformed under duality and conclude that the most relevant region (from the prac¬ 
tical viewpoint) is represented by the hrst quadrant of the (F B , F E ) plane. In section [5] 
the different portions of the (F B , F E ) plane are scrutinized in the light of the backreaction 
constraints (and for different initial conditions of the gauge couplings). The allowed region 
of the parameter space is explicitly obtained and discussed in detail. Section [6] contains our 
concluding remarks and a summary of the main findings. 


2 General parametrization of the coupling functions 


The general form of the coupling functions shall now be discussed hrst in the context of 
relativistic hydrodynamics and then in the case of conventional inhationary models. This 
analysis complements the parametrization already outlined in Eq. (1.4). 


2.1 Covariant decompositions 

To begin with let us recall the very well known covariant decomposition of V pu a namely: 

V fjUa vFUjj V ^yU a T (Ta/3 T T “T* a P’ (2-1) 

4 The duality transformations used here are a simple generalization of the strandard electromagnetic 
duality [23] . While in the standard case the electric and the magnetic gauge couplings coincide, in the 
present situation they differ. See section [2] for a more detailed discussion. 


4 




where 0 = V a u a and the projector is defined as V a p = (g a p ~ u a up). The remaining terms 
in Eq. (2.1) are given by: 

(Tap = ^Sp u a + V^) - ^ 7 [(Vyu a )t^ + ('Vr y up)u a \ - ^V a p, (2.2) 


UJ(x(3 2 ^P u OL ^ ’ otV*p) IX 7 [(V 7 ?Z a )?Z0 (V/y'Zi^)'Zi Q! ]. 


(2.3) 


In what follows, we shall assume that 0 J a p —y 0. Thus the symmetric coupling tensors A4 


pa 


and J\f P a can be written as: 


At pa A^t/pU^- + 

A“pa \bVj pUa T 


V f 


Y pU a 


2 M 

2mT”“‘ 


Vo-Mp ) + 


p^a 


VaUp 


V, 


2 M 

Z>b _ 7 

— -u 1 


-u ' 


(V 7 'Up)'U IT T (V 7 'U cr )'U / r 


2 M 


(V ryUp^U a T (Y jVia^Uf. 


(2.4) 

(2.5) 


where M denotes a generic mass scale; u p (and u p ) obey g pa u p Ua = 1 and g pa u p u a = 1- 
Notice that we dropped the terms containing g p „ since it can be reabsorbed by a redefinition 
of A in Eq. (1.1). 

In the analysis of [_22j the D-terms have been neglected, i.e. T>e = T>e = 0 and T>b = 
T>b = 0. All the terms appearing in Eqs. (2.4) and (2.5) will now be considered on equal 
footing. From Eqs. (1.1) and (1.2) the explicit form of the evolution equations of the gauge 
fields is: 

1_« 1. 


+ -Y a Z a ? - ^V q W q/3 = 4 tt/, 


( 2 . 6 ) 

Y a Y a/3 = 0. (2.7) 

Defining E a/3p< ^ = e a/3p< ’/ v / = ^ (where e a/3p ^ is the total antisymmetric pseudotensor of fourth 
rank) the two antisymmetric tensors Z a P and W“ /3 are given by: 

Z aP = M a ya(S _ M 0 y™ } = E^ Y a( Af p . (2.8) 

We shall now focus on the case where the background metric is conformally flat: 

gpu{r) = a 2 (r)rip V , 


n = ~, 

a 


(2.9) 


where r/ /t „ is the Minkowski metric and the prime denotes a derivation with respect to the 
conformal time coordinate r. Recalling that the standard Hubble rate is H = Ti/a, the 
explicit components of Z af3 and W al3 ar^J 


z^ = -z^ = -\\ e e^h 

r^B e 




V F e ljk 

Z lJ = -2 H——b k , 
M a 2 ’ 


M a 


V 


W h = — W iU = 2H= ( X B + ) ~—b k . 


M 


c ijk 


gt 


( 2 . 10 ) 

( 2 . 11 ) 


5 We recall that in the conformally flat metric of Eq. (2.9) the gauge field strength can be expressed in 
terms of the electric and magnetic fields as Y' ,a = e*/« 2 and Y l: > = —e l ^ k bk/a 2 . The explicit components 
of A ip and N p can be explicitly obtained from Eqs. (2.4) and (2.5) and they are: M® = A_e, A§ = A b, 
M\ = V E H6j/M and A? = V B H8 j JM. 
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2.2 Equations of motion and duality 


A eE + 4 nJ, 


Thanks to Eqs. (2.10) and (2.11), Eqs. (2.6) and ( |2.7 ) imply the following form of the 
equations of motion: 


( 2 . 12 ) 

(2.13) 

(2.14) 

(2.15) 

(2.16) 

(2.17) 


V x UA b BJ = S T 

v.f^Wo, 


V ■ (\/A 


where we have introduced the following rescaled variables: 

B = a 2 \J~X~b b, E = a 2 \[~X~e e, 
A - \ + Xb + ' Db h+ T>e h 

Ab - X+ ^2 + 2M H+ U H ' 

A - X 4. Xe + Ve H 4- Vb H 

Ae - X + T + 2M H + W H - 


= 47 TP, 


Concerning Eqs. (2.15), (2.16) and (2.17) three comments are in order: 


• in the homogenous case the terms containing T>e and T>b do not contribute to the 
evolution of the gauge fields; 

• in the simultaneous limit T>b —>• 0 and T>e —>• 0 we have that A b = A + A^/2 and that 
A^ = A + Ae/ 2, as previously established; 

• if we transform X E — > As, V B — >V E and V E —>> E> B then A^ —> A B (and vice-versa). 

The susceptibilities xe , and xb are defined as \e = x/Ae and xb = \/A b- The corresponding 
gauge couplings are instead the inverse of the susceptibilitie^] 


9 e — \/4vr/As, 9 b — \/4vr/As. 


(2.18) 


The mismatch between the two susceptibilities (i.e. / = Xe/Xb = ^-e/^-b) determines an 
effective refractive index affecting the evolution of the canonical fields. 

In Eqs. (2.12) and (2.14) J and p denote, respectively, the rescaled current density 


and the rescaled charge density. In what follows we shall analyze the amplification of the 
gauge fields from vacuum fluctuations so we shall set the sources to zero. There are however 
different possible initial conditions like the ones stipulating that initially a globally neutral 
plasma is present, namely p = 0 with J given by the Ohmic current of the plasma. This 


6 Note that in the limit Xb = Xe —> 0 and T>e = 'Db 
the gauge couplings, i.e. gE = 9b = \J 47t/A. 


0 xe and Xb coincide and the same happens for 
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possibility has been investigated in the conventional case ra bust shall not be specifically 
analyzed here. 

An important aspect that will be systematically used to simplify the discussion of the 


parameter space concerns the symmetries of the model. Neglecting the sources, Eqs. (2.12), 


(2.13) and (2.14) are invariant under the generalized duality transformation [22] that co¬ 


incides with the conventional duality m in the case of coincident electric and magnetic 
gauge couplings. When the gauge couplings are exchanged and inverted (i.e. g E —> 1 /g E 


and (]b —t 1 /g E ) the generalized duality transformation stipulates that Eqs. (2.12), (2.13) 
and (2.14) keep exactly the same form provided E —> —B and B —>■ E. 


2.3 Coupling functions and scalar fields 

We shall now examine form of the coupling functions when the four-velocity is related to the 
covariant gradients of a scalar field (i.e. either the inflaton or some other spectator field); in 
this case from the general considerations discussed above we shall have 


M pa (y) = p (fV a (p + + 


BM(. 

+ V 7 Vg-V? + Vg-t/7 V 7 V p y), 

KM) = + 


(2.19) 




V 7 ^ ptfj V 7 V cr -0 + V 7 V p -0^. (2.20) 

As discussed before M. pa and J\f pa do not have to depend on the same field. For illustration 


two different scalar fields (i.e <p and -0) have been used in Eqs. (2.19) and (2.20). The 


general form of Eqs. ( 2. 19| ) and (2.20) suggests that it is always possible to select the 
arbitrary coupling functions in such a way that the electric and magnetic gauge couplings 
depend on the scale factor. During the inflationary phase their analytic dependence can 
always be written as 


9e(p) = g E 

f(a) = fJ-Y, 

\a,iJ 


9s( a )=9B^'j > 


r 9 b __ 

B ~ -2 — ~T ■ 


9e 


where F E , F E and F denote, respectively, the normalized rates of variation, i.ej^] 


F e = 


din g E 
d In a 


F b = 


dlng B 

d In a 


F = 


d In / 
d In a ’ 


( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 


and, by definition, F = 2 (Fs — F E ). In Eqs. (2.21)—(2.22) Oj denotes the value of the scale 
factor at the onset of the evolution of the gauge couplings; this moment will be taken to 


'We recall that in this paper In denotes the natural logarithm while log denotes the common logarithm. 
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coincide, for simplicity, with the onset of the inflationary phascj^J Recalling Eq. (2.18), the 


parametrization of Eqs. (2.21) and (2.22) implies that x E (a) oc a Fe and \b oc 


1-Fb 


The parameter counting of Eqs. (2.21) and (2.22) goes as follows. There are nominally 


four parameters in Eqs. (2.21) and (2.22), namely the rates F E and F B and the corresponding 
amplitudes, i.e. g E and g B . These four parameters can be identified with F E , F B , and 
with the initial value of one of the gauge couplings (be it either g B or g E ). Neither g B or g E 
appear alone in the hnal electric and magnetic power spectra but only combined in f, which 
can be different from one and, in principle, also very large or very small depending on the 
specific model under consideration. On top of these three parameters there is also the total 
duration of the inflationary phase and, possibly, the overall duration of the pumping action 
of the gauge couplings. 


2.4 Power-law backgrounds 


Before proceeding further, the parametrization of Eqs. (2.21) and (2.22) shall now be justified 


in the context of specific scenarios. This step is, in a sense, superfluous given the arbitrariness 
of the coupling functions. However it might be interesting to have this point treated in more 
depth. In what follows we shall assume the validity of the standard inflationary evolution 
and examine different specific models. The power law backgrounds are interesting insofar as 
they are exactly solvable. In this case it is immediately evident that the coupling functions 
must have an exponential or quasi-exponential form if we want to justify the validity of the 
parametrization of Eqs. (2.21) and (2.22). In power-law backgrounds e = — H/F [ 2 = 1/a 
and 


9 2aAF p 


V = <Po + \ ~ M P In (t/U), a(t ) = (t/ti) c 


(2.24) 


where t t denotes the initial cosmic time of the inflationary expansion. It can also happen 
that the susceptibilities do not depend on the inflaton but only on some spectator field that 
evolves during a quasi-de Sitter stage of expansion. In this case the evolution equation of ip 
can be solved even exactly like in the case of an exponential potential. If the field ip evolves 
in an exact power-law background its evolution will be given as 


ip = ipo + F X M P In (t/U), W(ip) = Wo e- K2 ^-^ Mp , 


(2.25) 


_2 

where /C 1 /C 2 = 2 and R 0 /C 2 = H/M P e(e — 3). With this information the explicit couplings 


can be explicitly constructed and they will fall within the parametrization of Eqs. (2.21) 


and (2.22) 


I 11 the conventional case the coupling functions are exponentials mmm - if we make 
the same choice here A(y>) = Qo exp [qo(g> — </?o)] sets the overall normalization while the 


8 If this identification is not made there will be a supplementary parameter fixing the origin of the evolution 
of the gauge couplings in comparison with the onset of the inflationary phase. 

















other couplings can be parametrized as: 


Bi(<p) = Q 1 


*qi(<p-<po)/Mp 

HMl ' 


B 2 {<p) = Q 2 - 


<Po)/M p 

H}Mp 


BM = Q 3 - 


^q 3 (<fi-tpo)/M P 


HfMr 


(2.26) 

Exactly the same parametrization can be envisaged for C\(ip), C 2 (ip) and C 3 (ip) where we 
choose, for simplicity, ip = p: 


e Pi(<p-tpo)/Mp 

C\Up) = Pi - —o —. 

II? Ml 


e Pa(<p—<po)/Mp 

CM = p 2 -==—. 

HfMp 


e P3(<p—<Po)/Mp 

c 3 Ip) = P's —p=3—■ 
HfMl 


(2.27) 


Equations (2.26) and (2.27) extend the class of models of Refs. pnnununi to the case of 
derivative interactions obtained through the present construction. 

The explicit form of the couplings of Eqs. (2.26) and (2.27) together with the background 
evolution of Eq. (2.24) leads to the parametrization of Eqs. (2.21) and (2.22). While 
more complicated dynamical situations are certainly possible (and should be separately 
investigated) we shall now argue that the monotonic evolution of the gauge couplings is 
realized in various explicit scenarios. Let us assume, as above, ip = p and also set = 

B 2 (p) = B 3 (p) = 0 together with C 3 (p) = 0. The explicit form of the susceptibilities, in the 
case of exponential couplings, will then 


A i?0) 

— Qo 




V r*/ 

A sO) 

— Qo 



V r*/ 


P 2 
+ w 

v^e(l 


P 2 V2e( — 


T \ 


V2e(p2-V2e)/(e—l) 


T, 




(2.28) 


r 


\ \/ 2 e(p 2 —\/2e)/(e—1) 


T \ C 2 e(p!-\/ 2 e)/(e-l) 


(2.29) 


Depending on the relative hierarchy of q 0 , pi and p 2 the paranetrization of Eqs. (2.21) and 
(2.22) can be easily recovered after some putative time related to r*. This parametrization 
may also lead to more complicated situations where /(r) may even have a non-monotonic 
behaviour. However it is sufficient to suppose that p 2 q 3 and p\ 3> p 2 qo to obtain 


XE(a) = Xe(J^) , Xb{o) = 


-f b 


(2.30) 


where F E = \Je/2{p 2 — \/2e)/(e — 1) and F B = \Je/2(pi — y / 2e)/(e — 1); moreover xp = \/2 eP 2 
and xl = e Pi- The initial value of the susceptibilities and the relative strength of the 
couplings can be tuned by changing P 3 and P 2 . 

9 The connection between cosmic and conformal time can always be written as (t/tp) = (—t/t*) 1 ^ 1 -") 
where r* = U/(a — 1). 
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2.5 Slow-roll models 


The previous discussion can be easily generalized to the case of generic slow-roll evolution. 
The simplest situation in this respect is given by the case where all the undetermined func¬ 
tions vanish except two of them. Let us therefore set Bi(<p) = B 2 (<p) = B 3 (p) = 0 and also 
C 3 (<p) = Ci(<p) = 0. Given a certain inflationary potential we can always define 


Z(<Pi,<p) 



W 

w- 


dp. 


The coupling functions can then be fixed as: 


C 2 (p) 

\(ip) 


6 


e -2 F E X(tpi,ip) 


g— 2 (F B -F E )l(tpi,tp) _ 


g-2 F E Z(tpi,y) 


2g —2( F B—FE)Z(<Pi,<P) _ 


Another possibility is to set C 2 (<p) = C 3 (p) = 0 and, as before B\(p) = B 2 (p) 
The values of A(</?) and Ci(p) are given, in this case, by: 


(2.31) 


(2.32) 
B 3 {p) = 0. 


X(p) = Ci(^) = -— e - 2 W<Pi.v) 


g—2 (FB—F E )T((pi,p) _ 


(2.33) 


Equations ( |2.32 ) and (2.33) imply that the susceptibilities can be written, in a unified no¬ 
tation, as: 


Xb 


= A + hMy_ C2 („) 


H<p+ ^ 


= Xb 


X 2 e = A + C 2 (<p)Htp = x i E (j^j 


— 2 Fb 


—2 Fe 


(2.34) 


where we recall that in the slow-roll approximation (a/a*) = exp [—I(</?j, </?)]; furthermore 
the slow-roll relations have been used where appropriate. The same analysis leading to Eq. 
(2.34) can be repeated when all the Ci(ip) are equal to zero while, on top of A(</?), the only 
two non-vanishing coupling functions are either Bi(ip) or B 2 (p). 

It is obvious that depending on the model the relation between (a/aj) and tp will be a bit 
different but still the general parametrization of Eqs. (2.21) and (2.22) will remain valid. In 
the case of simple monomial potentials we have W = and therefore we shall have that 


= exp 




2 pM P J 


(2.35) 


Other typical models can be analyzed like the case of small field and hybrid models where 
the potential can be written as W(p) = M 4 ( 1 ± tzp v ) (where the plus corresponds to the 
hybrid models while the minus to the small field models). In the case of R 2 models the 
potential, in the Einstein frame, is given by 

H-'W = 3M ^ ( 1-e-^/M,) 2 , (2.36) 
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where we defined, for practical reasons, = \J2/3tp/M P . Using the evolution equations we 
obtain: 

$ = In [2M(f* — f)/3], a = a* exp 3[<f> — e $ ]/4. (2.37) 

In cosmic time the scale factor reads 


a{t) = [M(U-t)} 3/4 e~ M ^- t)/2 . 


(2.38) 


where t* denotes in practice the final time of the inflationary expansion. All the specific 
analyses discussed before can be easily applied also to this case with just one caveat. To 
obtain the specific equations relating the cosmic and the conformal time parametrization we 
must expand the above equations in the limit i C 1,. The result of this simple expansion 
can be written as 


a~a*x 3 J 4 e- x */ 2 eP x - 3 W x *W 4 : 


(2.39) 


where x — Mt and = Mt *. Thanks to the previous expression we can easily deduce the 
relation of t to r by recalling that a(r)dr = dt. In summary we have demonstrated that 
while the most general parametrization of the coupling functions may lead to non-monotonic 
susceptibilities, after some initial transient time the monotonic behaviour can be recovered. 
The parametrization of Eqs. (2.21) and (2.22) is therefore fully justified and it encompasses 
a variety of models. 


3 Quantization and power spectra 

The problem is greatly simplified by abandoning the conformal time and by rescaling the 
canonical variables. This transformation leaves the action invariant and simplifies the evo¬ 
lution of the mode functions. 


3.1 Time reparametrization of the action and quantization 

The main idea is to use a new time parametrization in the action and a new set of rescaled 
vector potentials. In time-dependent (conformally flat) backgrounds and in the Coulomb 
gauge (i.e. Y 0 = 0 and V • Y — 0) that is preserved (unlike the Lorentz gauge condition) 
under a conformal rescaling of the metric the action is 


S = 


dr d 6 x \ A' 2 + 


XeY 

Xe) 


/ 2 

A 2 - 2■ A' - diA ■ d { A k 

Xe Xe 


(3.1) 


where 10 


A = JA.e/(^e)Y. We have assumed that Xe and \b are only dependent on the 


conformal time coordinate r. In terms of the canonical momentum conjugate to A the 


10 The l/VYr is purely conventional and its presence comes from the factor 167 t included in the initial 
gauge action. 
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canonical Hamiltonian is simply given by [22]: 


Ha(t) = - J d 3 x 


x'e -> r &iA ■ d l A 

n z + 2^7f • A + ——-- 

Xe f 


j? = A' - — A 
Xe 


We can now change the time parametrization in the action of Eq. (3.1) 


T ^r] = dr = 

and simultaneously redefine the vector potential: 

A(t,x) A( V ,x) = JTJ^yy 

Ultimately A can be directly expressed in terms of Y as: 

A E {rj) A B (rj) 


(3.2) 


(3.3) 


(3.4) 


A(rj, x) = 


Ae(v) 


\ ^\/7W) 


Y(v,x) = 


\[Yk 


Y(r],x) = 


Y(j],x) 


9e(v)9b(v) 


(3.5) 


where the second and third equalities follow, respectively, from the mutual relations among 
the different susceptibilities and from the relation between each susceptibility and the cor¬ 


responding gauge coupling. Using Eqs. (3.3) and (3.4) the action (3.1) becomes: 


S=\f\[ dr] d 3 x\ a ? X d n A+ f ^ XeXb L \ 2 J2 _ 2 A-A- &A- B 1 a\. (3.6) 

2J 2J l L y/X eXb \ x/XAXb J 

where the thick overdot (or bullet) denotes a derivation with respect to rj. From the action 
(3.6) the new canonical Hamiltonian and the corresponding canonical momenta are now 


given by: 


h a(v) = 2 J d * x 

fl(r],x) = d v A — 


n 2 + 2 


(VXeXbY 

\JXeXb 
(VXeXbY r 


■U-A + diA-tfA 


\JXeXb 

The Fourier mode expansion for the canonical fields: 
_ I _ [ a3 

(2tt) 3 / 2 


(3.7) 

(3.8) 


A(x,v) = fr>S\Z/9. I d 3k HfcCh) e 


—ik-x 


A&V) = j k Mv) e 


—ik-x 


(3.9) 


implies that the canonical Hamiltonian (3.7) can be rewritten as: 
1 


Ha(v) = 2 Jd 3 k 


n £ • n_ E + { -^= (U • A-t + n. E • A c )+eA- k -A_ i 


(3.10) 
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From Eq. (3.10) the corresponding equations or motion are: 


+ \Z9b9e 


V9b9e 1 k ’ 


drAk ~ ~h 2 Aj: — ^9b9e 


n r 


(3.11) 

(3.12) 


\Z9b9e- k ’ 

where we have used the explicit relations connecting the susceptibilities to the gauge cou¬ 
plings. The duality transformation exchanges the canonical fields and the conjugate mo¬ 


menta so that Eqs. (3.11) and (3.12) go one into the other by virtue of the following 
transformation: 


Mv) -* 3bM -* aka' 


n t M -> -kA s ( v ), Mv) -> (3.i3) 


Promoting the canonical fields to quantum operators (i.e. A t —> Ai and 7 q —> n,) the 
following (equal time) commutation relations (in units fi — c— 1 ) must hold: 


[Ai(x u 7}rj)} = - x 2 ), A„(fi - x 2 ) = 


d 3 k 

(2vrj 3 


(k), (3.14) 


where Pij(k ) = (8 l3 — kAj/k 2 ). The function A tJ (x\ — x 2 ) is the transverse generalization of 
the Dirac delta function ensuring that both E and A are divergenceless. The held operators 
can then be expanded in terms of the corresponding mode functions 

d 3 k 


Ai(r), x) = 


T.A(k) 


(27t) 3 / 2 „ 

d :i k 


F t (v) * + Fl(rf) al „e iki 

t 


n*(77, x) = J 2 51 e i\k) G k {r)) d Ka e lk ' x + G* k {rj) a£_ a e' 
where F k {rj) and G k (r] ) are the mode functions in the 77 -parametrization. 


(3.15) 

(3.16) 


3.2 General solutions for the mode functions 


As a consequence of the indetermination relations, Eqs. (3.15) and (3.16) imply that the 


mode functions must obey the equations derived previously and the Wronskian normalization 
condition: 

F k ( v )Gl( V )-Fl(v)G k ( V )=z. (3.17) 


From Eqs.(3.11) and (3.12), F k (rf) and G k {rj) can be shown to satisfy the following pair of 
equations: 


dF i, ^ 


drj 

dG k 

drj 


— G k + a J~9b9e[^ 


1 V F k , 


\/9b9e 


— ~k~F k — sj g b9e 


a/ 9b9e 


G 


k- 


(3.18) 

(3.19) 
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Combining Eqs. (]3.18|) and (|3.19|) we obtain two decoupled equations: 

1 


d 2 F k 
drf 
d 2 G k 
drf 


+ 


+ 


k — GQbQe 


\/9b9e 


F k = 0, 


k 2 - 


(a J9b9e)' 


\/9b9e 


G k — 0. 


(3.20) 

(3.21) 


Equations (3.20) and (3.21) can also be written in terms of the susceptibilities xe and \b by 
simply taking into account that, as already mentioned, Qe = \JFeJxe and (]b = \J^/xb- 
Equations (3.18)—(3.19) and (3.20)—(3.21) can be directly derived from the mode functions 


obtainable from the quantization in conformal time. Indeed, using the recipe 
dr 


v / 7FM 


= dri 


Fk(r) = / 1/4 [^(T)]F fc [r/(T)], G k [ri(r)] = / 1/4 [h(^)]G fc [?/(r)], 

(3.22) 


Eqs. (3.18)—(3.19) become immediately 


Z 7 >/ _ r 1 A ~^-B 77 

r u — Lt£. I - -^k') 

Xe 


G' k = 


e 

f 


Xe 

Xe 


G k , 


(3.23) 


where the prime denotes a derivation with respect to r; Eqs. (3.23) coincide with the ones 
derived in 


The general solutions of Eqs. (3.20) and (3.21) can be obtained either exactly (i.e. by 
solving the evolution in exact terms) or approximately (i.e. in different asymptotic limits). 


In a large class of physical situations the pump holds of Eqs. (3.20) and (3.21) are expected 


to decay for \rj\ —> oo faster than l/r] (see however at the end of this subsection for a different 


situation). The solution of Eq. (3.20) can then be obtained as: 
1 


F k(v) = 


zb ikrj 


V2k 


F k (v) = 


D 


(i) 


+ 


D 


k » 


( 2 ) 


V 9 e9b 


\Z9b9e 


\Z~9e~9b \[9e 9 b 


dVi9E(Vi)9B(Vi), 


k 2 < 


\Z9e9b 


V 9 b9e 


(3.24) 

(.3.25) 


where D ^ and D^ are the two arbitrary constants that must be fixed from the boundary 


conditions. Recalling Eq. (3.18), Eqs. (3.24)—(3.25) also determine the explicit form of 


G k (r]). Alternatively Eq. (3.21) can be directly solved with the same strategy leading to 


Eqs. (3.24) and (3.25). 


3.3 Exact solutions with asymmetric gauge couplings 


According to Eqs. (3.20)—(3.21) and using the generalized parametrization of Eqs. (2.21) 


and (2.22), the evolution of the mode functions obeys: 
a 2 - 1/41 


d 2 F k 
drf 


+ 


k 2 


V 


F k — 0, 


d 2 G k 

drj 2 


+ 


k 2 



G k — 0, 


(3.26) 
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whose solution is is given by: 


F ^) = 4^ 


1-2 Fp 


V2k 


a = 


2(1 + Fb — Fe) 


Ik 


G k (r,) = -M\llsl-kr ] H^l 1 (-kri), 


(3.27) 

(3.28) 


where |A/”| = \Jtt/2. The case Fe = Fb + 1 is a bit special since we have that / = /(r/rj) 2 . 
But in this situation the susceptibilities can be parametrized, in general terms, 

Xe( if) = XB (r,) = Xge «B(i+ml/i % /(,) = J e 2(o«-o.)(i+n)/i.. ( 3 . 29 ) 


The evolution equation for F k {rf) can be written as: 


d 2 F k 

dr] 2 


+ 


k 2 


(Ol e + OIb ) 2 


F k 


0, 


and the solution of the previous equation is: 


Fktn) 


1 

Vxh e 


iQ k (r]+rji) 

9 





(ats + as) 


4r/ 2 


(3.30) 


(3.31) 


If kr /* 1 then Q k ~ k while in the opposite limit £2^ ~ z|oe + ob|/( 277 *). Recalling the 

expression of G k we have that 


Gfc — —z 



*(ofE + Ofs) 
2 £ 277 * 


C+v*) 


(3.32) 


Note that in the limit fcr/* < 1 we have that G k {j]) —> 0 and PE(k,r)) —> 0 provided and 
Ob are both positive semidefinite. Conversely, in the same limit, 


n(^r 


4?y* 

ct£ + as 


g|os+ Q s I ( r H _7 7i)/ J 7* 


(3.33) 


3.4 Correlation functions 


In the r/-pararnetrization the correlation functions of the electric and of the magnetic fields 
can be swiftly computed. It is sufficient to recall that the canonical fields can be directly 
expressed in terms of the comoving electric and magnetic fields, i.e. E and B: 


E(v, x) = 


gfaj) 


B(r],x) 


V x 


-A(ri, x) 


(3.34) 


n In this specific case the relation between r and r/ is exponential since rj is the logarithm of r, i.e. 
In (-r/rd = - V/ffi/ffi) and % = Ti/sf]. 
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In terms of the mode decomposition the comoving electric and magnetic fields are 


Bi(v,x) = 


ie r 


Y / d 3 kk m e ( Y F k (g)d?e 


—ik'X 


t 




(27t) 3 / 2 a 

EA r/,x) = - 1 ._ Y [ d 3 ke^\G k (ri)dre- i]; - s + Gl(ri)at e*' 

V ’ (2vr) 3 / 2 V V ’ k ’ a ky ’ fc ’“ 


,(3.35) 
. (3.36) 


In terms of the mode functions, the Fourier components of the field operators Bi(x,r) and 
Ei(x,r ) are respectively: 


Q) — r Cmni ^ Qm[dq,a Fq(v) + E q (r])\, (3.37) 

y/fW 

Ei(v, q) = -7= Y e i a) [dq ,0 G q {rj) + d)_^ G* q (rjj). (3.38) 

tffiv) “ 

The correlators in Fourier space are then given by: 

- 27T 2 

{Bi(r),k) Bj(ri,p)) = -p- P B (k, rj) P^k) 5 {3 \k + p), (3.39) 

- 27T 2 

{Ei(r), k) Ej(r),p)) = — P E (k,rj) P^k) 5^ 3 \k + p), (3.40) 

where Pij(k ) = (S,j — kjkj/k 2 ) and where the power spectra are: 


PB{k,rj) = 

k 5 

\F k (v )\ 2 , 

(3.41) 

2tt 2 a 4 ( j]) \Jf{j}) 

PE(k,v) = 

k 3 

\G k (v)\ 2 - 

(3.42) 

2vr 2 a 4 {rj)Jf{ri) 


Note that, as a consequence of the duality symmetry, the magnetic and the electric power 
spectra are interchanged when g E — > 1/gB and gs — > 1/gE • In this case G k —> —kF k and 
F k —> Gk/k so that Ps(k, 7 ]) PE(k,rj ) and vice-versa. It is very important to remark, for 
the following discussion, that the definition of power spectrum of Eqs. (3.39) and (3.40) is 
particularly suitable since P E (k, rj) and P E (k,r]) have both dimensions of an energy density. 
Consequently we shall measure the amplitude of the power spectrum either in critical units of 
in nG 2 . We stress that this occurrence is not only mathematically useful but also physically 
sound: if the power spectrum is defined according to the present strategy it is evident, as it 
has to be, that the magnetic power spectrum is also (up to an irrelevant numerical factor) 
the magnetic energy density per logarithmic interval of wavenumbers (or logarithmic interval 
of frequency). 
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4 Power spectra and initial data 


4.1 General form of the power spectra 

The power spectra depend on the initial data of the electric and magnetic susceptibilities 
(or, which is the same, of the electric and magnetic gauge couplings). Two complementary 
possibilities can be envisaged. The first possibility stipulates that /(a*) — fi — 1, he. 
Xs(oj) = Xs(ai): the susceptibilities are initially equal and then they diverge so that, 
depending on the number of efolds, they will still be different at the end of inflation. The 
second possibility is the one where /(a/) — ff — 1; in this second case xb(o/) = Xs(a/), he. 
the susceptibilities are initially very different but they become equal at the end of inflation. 
In this second case /(a*) 1. The various intermediate situations (stipulating that initially 

there is an arbitrary mismatch between the susceptibilities) will fall between the patterns of 
the two limiting cases f, — 1 and ff — 1. 


With this distinction in mind, and recalling Eqs. (2.21)—(2.22) and (2.23), the solution 


M-i ( k \ 
! \aH) 


P B (k , a, cr) = H 4 Q b {<j,h) /, 

Qb(<t,//) = 1—hhliZ 2 2 M -3 |i -)_ /Lt | 2 l cr l— 1 


k \ 5—2|ct| 


of the mode function of Eqs. (J3.27|) implies that the magnetic power spectrum becomes: 

a \ 2/4M- 1 ) 

r 2 (M) 


(4.1) 


7T 0 


where, for practical reasons, the following notation has been employed: 


1-2 Fp 


a = 


1 + F b — Fj 


H ^ — F B — F E . 


Similarly thanks to Eq. (3.28) the electric power spectrum is: 

/ U \ 5—2\a — 1| 

P E (k, a, cr) = H x Q E ((J,/i) /j CT " 1|_1 J 


2/xdcr—1|-1) 


(4.2) 


(4.3) 


Qe{ct, h) = 


r 2 (l^-i|) 


7T 


2 2k-H-3| 1+ ^|2k-l|-l_ 


Note that in the plane (F B , F E ) there is a singular trajectory, namely 1 + F B — F E = 0 where 
cr diverges. This singularity is not physical and stems from the fact that for F E — F B + 1 
the gauge couplings evolve exponentially in r/. 


An essential feature of Eqs. (3.39) and (3.40) is that the electric and the magnetic power 


spectra both have dimensions of an energy density in Fourier space and this is related to the 
fact that the scale-invariant limit of the spectra occurs for n B —y 1 and n E —> 1 . As already 
mentioned we shall measure the magnetic power spectrum either in critical units of in nG 2 . 
Another way of measuring the amplitudes of the power spectra will be to use critical units 
and, in this case, the magnetic and electric power spectra can be expressed as: 

87 t 2 ii-, / k \ 5_2 l cr l / a \ Mkh 1 ) 

n B (k, a, a, /i, ft) = — e AtzQb(ct, /i) /j ff| J J , (4.4) 
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Qn f2 / U \ 5— 2\a— 1| 

Q E (k, a, a, /i, fi ) = — eUl^ /i) /f“ “ 


(3-2 


2/x(|o—1| —1) 


. (4.5) 


Notice that in Eqs. (4.4) and (4.5) we traded H/Mp for A-p, i.e. the amplitude of the power 


spectrum of curvature perturbations at the pivot scale k p = 0.002 Mpc ; note also that the 
slow-roll parameter e appears in the final expressions. 


The magnetic and the electric spectral indices can be easily deduced from Eqs. (4.1) (4.3) 


and from Eqs. (4.4)—(4.5) and they are: 


n E = 6 — 21 cr |, n E = 6 — 2 |cr 


(4.6) 


If a certain hierarchy between the gauge couplings is imposed ab initio for some theoretical 
consideration, then J is fixed. Similarly fi can also be fixed by tuning /(a/) at a certain 
specific value, be it ff: in this case ff = f/af/af) 2 ^. In the conformal time parametrization 
the spatial gradient leads to a term k/y/J which can be either blueshifted or redshifted 
depending on y/J . Naively we would require that y/J > 1 to avoid a light speed that is larger 
than one. This observation, however, depends on the time parametrization. For instance in 
the cosmic time parametrization we have that the physical frequency is given by k/a(t ); if 
we would apply the same tenets we would exclude the interesting possibility of contracting 
Universes. In the specific cases discussed below the cases y/J <C 1 are phenomenologically 
disfavoured and are anyway excluded for more physical reasons. Let us finally note that to fix 
fi or ff is in general different. However in the framework of the monotonic parametrization 


of Eqs. (2.21) and (2.22) the value of ff is uniquely determined as a function to J through 
the total number of inflationary efolds, i.e. In (////*) = 2 (Fp — F E ) N t . Therefore, for the 
purposes of the concrete discussion, the boundary conditions can be fixed either in terms of 
fi or in terms of ff. 


4.2 Physical regions of the parameter space 

The 6 different regions of the parameter space in the (F E , Fp) plane are illustrated in Fig. [I] 
where on the horizontal and vertical axes we report, respectively Fp and F E , i.e. the rates of 
variation of g B and g E . The different regions appearing in Fig. ([Tj) correspond to 6 different 
dynamical situations that can be summarized for short as follows: 

(I) in the first region the gauge couplings are both increasing, i.e. Fp > 0 and F E > 0; 
furthermore also / increases since Fp > F E ; 

(II) in the second region the gauge couplings are also increasing but since F E > Fp we 
have that F < 0 and / is bound to decrease; 

(III) in the third region Fp < 0 and F E > 0 (leading to F < 0); 

(IV) in the fourth region both gauge couplings are decreasing, (i.e. Fp < 0 and F E < 0) 
while Fp < F e and consequently F < 0; 
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Figure 1: The parameter space of the model is schematically illustrated in terms of the 
normalized rates of variation of the gauge couplings. 


(V) the gauge couplings are also decreasing simultaneously (i.e. F E < 0 and F B < 0) but 
in this region F > 0 (i.e. / increases); 

(VI) in the sixth and last region F E < 0 and F E > 0 implying F > 0 and the corresponding 
growth of /. 

In regions (I) and (II) the gauge couplings are both increasing; in regions (IV) and (V) 
the gauge couplings are both decreasing (but at a different rate); finally in regions (III) and 
( VI) one of the two gauge couplings is increasing while the other is decreasing. Recalling 
Eq. (4.2) the spectral indices can be directly expressed in terms of F E and F B : 

1-2 F e , „ „ x „ 1 + 2F, 


n B (F B , F e ) = 6 — 


(1 + F b — F e ) 


n E (F B , F e ) — 6 — 


B 


(F e — F b — 1) 


Thanks to Eq. (4.7) it is immediate to show that, under duality: 

F e —> —F b , F b —» —F e , n B —> n E , n E —> n B . 


(4.7) 


(4.8) 


Furthermore, recalling Eq. (4.2), under duality a —>• a = 1 — a. 


The magnetic power spectra solely expressed in terms of F B and F E become then: 

r . , „ „ m /n / k \n B (F B ,F E )-l 

P B (k, F B , F E , fi) = H 4 Q B (F B , F E )ft nB(FB ' FE)]/2 J Af B (N t , F B , F E ), 

= r"[(6 - n B (F B ,F E ))/ 2] 2 3 -n B (F B ,F E ) M _|_ _p F \5-n B (F B ,F E ) 


7T° 

— e Nt{F B —F E )[i—n B {F B ,F E )] 


Q b (f b , f e , 

Nb(F b , F e ) 
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As a function of F E and Fe the electric power spectra become instead: 


P E (k, N„ F b , F e ) = H‘ Q e (F b , F, 


)fl i - n -^W(F) ns ' Fa ' FEyi M E (N t , F b , F e ) 


Qe(F b , F e ) 
A f E (N t , F b , F e ) 


\aH J 


r [(6 Ue(Fb, Fe))/2] 23 -n E (F B ,F E ) _|_ p^ _ p^ 5 -n E (F B ,F E ) 


7T 

= g N t(F B —F E )[4—n E (F B ,F E )] 


(4.10) 


where n E (F E , F e ) and n E {FB,F E ) are now given in Eq. (4.7). With the same strategy the 


corresponding critical fractions of Eqs. (4.4) and (4.5) can be easily expressed in the (F E , F E ) 


plane. Let us finally mention that Eqs. (4.9) and (4.10) go one into the other thanks to a 


duality transformation. This important property, already mentioned in the previous section, 


is now explicitly verihable at the level of the final results and thanks to Eq. (4.8). 


4.3 Dualizing the regions of the parameter space 

In each region of Fig. [lj depending on the values of F E , F E and F, the power spectra have 
different slopes and different amplitude.This aspect is illustrated in Fig. [2] where we analyze 
the parameter space from the viewpoint of the slope of the power spectra. The shaded 



Figure 2: Dependence of the power spectra on the rates of the gauge couplings. 

regions of Fig. [2] correspond to different spectral indices and different values of a. Under 
duality the various regions are interchanged. Consider, for instance, the first quadrant of 
Fig. [2]where the gauge couplings are both increasing, i.e. F E > 0 and F B > 0. In this region 
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Normalized rates 

Ur 

n E 

0 < F e < 1/2 

n-s = (5 + 6 Fb — 4Fk)/(l + Fb — Fe) 

he — 10 — ns 

1/2 < Fe < Fb + 1 

ns = (7 + 6 Fb — 8Fe)/(1 + Fb — Fe) 

n E = n B - 2 

Fe > Fb + 1 

ub — ( 5 + 6 Fb — 4:Fe)/(1 + Fb — Fe) 

he = n B + 2 


Table 1: First quadrant Fe > 0 and Fb> 0. 


we have: 


0 < F E < 1/2, 0 < a < 1, 

1/2 < Fe < Fb + 1, cr < 0, 

Fe > Fb + 1, cr > 1. (4.11) 


For each region of Eq. (4.11) the spectral indices have a different functional dependence 
upon Fe and Fb- In Tab. [l] the spectral indices in the first quadrant are summarized. 
Notice that the electric spectral index is always expressible in terms of ns but the relation 
changes from region to region. In Fig. [2]there are peculiar lines dividing the different shaded 
areas. The straight line F E — F B + 1 defines a line of singular models separately discussed 
in Eqs. ( |3.29[ )-( 3.30[ ) . These are not physical singularities but rather singularities of the 
parametrizatioij^j 

Always in Fig. [2] the lines Fe = 1/2 and Fb = —1/2 are the zeroes of a and of 1 — a. 
Both models are regular, the corresponding spectra have logarithmic corrections either in 
the electric or in the magnetic case. For the sake of simplicity, we shall exclude these three 
lines when charting the area of the parameter space in the plane (Fb, Fe)'- this is not crucial 
since we are just excising three lines from a whole plane. Furthermore these cases have been 
separately discussed in the previous sections, they are not more significant than the others 
and they are also irrelevant for the phenomenological considerations discussed in section [5] 

In Fig. [2] there are finally two particularly simple cases, namely F E = 0 (with Fb free to 
vary) and Fb = 0 (with Fe free to vary). These two situations correspond to the situation 
where one of the gauge couplings is constant. 

Different portions of Figs. [I] and [ 2 ] can be easily related via duality. Equation (4.11) 
refers to the first quadrant of Figs. [T| and [2} Let us now suppose to be interested in the 
electric and magnetic spectral indices in the third quadrant where F E < 0 and Fb <0. In 
this region both gauge couplings are decreasing. The three regions of Eq. (4.11) transform 


12 As previously shown in this case the susceptibilities have an exponential dependence on 77 and, con¬ 
sequently, the relation between r and 77 is logarithmic. These models are characterized by an electric field 
than vanishes exactly but besides this they do not have a particularly relevant status from the viewpoint of 
phenomenology. 
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Normalized rates 

riR 

n E 

0 < F b < 1/2 

Us = 10 — 77 -e 

tie — (5 — 6 Fe + 4Fg)/(l + Fb — Fe) 

1/2 < Fe < Fb + 1 

n B = 2 + n E 

n e = (7 + 8 Fb — 6Fe)/(1 + Fb — Fe) 

Fe > Fb + 1 

n B = 2 7i E 

tie — (5 — 6 Fe + 4Fg)/(l + Fb — Fe) 


Tabic 2: Third quadrant Fe < 0 and Fb < 0. 


under duality into: 


- 1/2 < F b < 0 , 0 < 

— 1 + Fe < Fb < — 1 / 2 , 

Fe > Fb + 1 , cr < 0 . 


a < 1 , 
a > 1 , 


(4.12) 


The result of Eq. (4.12) can be easily derived by recalling that, under duality, Fe —>■ —Fb 
and Fb —$■ —Fe- Notice also that duality relates the first and third quadrant (region by 
region) given the transformation properties of a. The results of the duality transformation 
are reported in Tab. [2j In the remaining two quadrants the same logic applies so that the 
spectra, the spectral indices and the other information can be easily deduced by successive 
duality transformations from the first or from the third quadrant^} 

Before discussing in detail the parameter space of the model let us remark that there are, 
broadly speaking, three distinct physical situations. The gauge couplings may be both small 
at Oj and then increase: this is the case of weakly coupled initial conditions evolving towards 
strong coupling at the end of inflation. The second possibility is that the gauge couplings 
are initially of order 1 and then decrease towards the end of inflation: this is the case of the 
strongly coupled initial conditions. The third possibility is that initially one of the gauge 
couplings is strong while the other is weak: this is the case of mixed initial conditions. In 
the case of conventional inflationary models the Universe evolves from strong gravitational 
coupling to weak gravitational coupling, i.e. the space-time curvature is maximal at the 
onset of inflation and gets smaller during reheating. It is fair to say that the potential 
drawbacks of magnetogenesis coincide with the potential drawbacks of conventional models 
of inflation which are, typically, not geodesically complete in their past history and preceded 
by a high curvature regime. The considerations reported here can be easily extended to the 
case of bouncing models evolving from weak gravitational coupling to strong gravitational 
coupling, i.e. the space-time curvature is small initially and gets larger at the reheating. 

In summary, it would be pedantic to discuss in details all the 6 regions of the parameter 
space of Fig. [I] (or Fig. [2]). So it is mandatory to complement Figs. [I] and [2] with more 
physical considerations. Fortunately duality comes to our rescue since the different portions 
of the parameter space are all related by duality transformations so that, for a complete 


13 Notice that from Tab. [I] the structure of Tab. [ 2 ] can be immediately obtained by only using duality. 
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discussion of the parameter space, a thorough examination of a single quadrant will suffice. 
Since the possibility of having the simultaneous growth of the gauge couplings during inflation 
is typical of this type of model we shall therefore focus on the first quadrant where F B > 0 
and F b > 0. 


5 Charting the parameter space of magnetogenesis 

The parameter space of inflationary magnetogenesis will now be charted by discussing in 
detail the first quadrant of Figs, [l] and [2j Thanks to the duality symmetry the remaining 
portions of the parameter space will be easy to discuss (see section [4] for a specific discussion 
of this point). Three independent requirements are relevant for the present discussion: 
the critical density bound (sometimes dubbed backreaction constraint), the magnetogenesis 
requirements and the naturalness of the initial conditions. All the mentioned conditions are 
necessary but none of them is per se sufficient to pin down the allowed region of the (F B , F E ) 
plane. However, if they are all jointly verified in a given domain they become sufficient. 

The logic will be, in short, the following. We shall first identify and classify all the 
domains of the (F B , Fe) plane where the backreaction constraints are satisfied. Then, in 
the regions not constrained by backreaction we shall compute precisely the magnetic power 
spectrum by keeping track of all terms and prefactors of the spectra. It will then be easier to 
see in the (Fe, F b ) plane where and how the weakly coupled initial conditions are verified. 


5.1 Maximal wavenumber of the spectra 

We recall that kr can be expressed as: 


kr 


k 

(1 -e)aH 


k 

Ho 


e 


N max 


l + e + 0 (e 2 ) , 


(5.1) 


where Hq = 2.334 x 10 -4 (ho/0.7)Mpc _1 is the present value of the Hubble rate and lV max 
is the maximal number of efolds which are today accessible to our observations [23] (see also 
|26j). In practice IV max is determined by fitting the redshifted inflationary event horizon 
inside the present Hubble radius H ( J ] : 


e v max = ( 27re ^p R0 ) 1 /4 


M P 

Ho 


V 2 (H r y-v* 

h) 


(5.2) 


where f 2 p 0 is the present critical fraction of radiation (in the concordance model — 

4.15 x 1CT 5 ). The term containing 7 accounts for the possibility of a delayed reheating 
ending at a putative scale H r eventually much smaller than the inflationary curvature scale 
F[. For illustration we shall focus on the simplest case and choose the sudden reheating 
approximation by setting 7 = 1 / 2 . Different possibilities can be however considered. A 
long postinflationary phase dominated by a stiff equation of state has been examined in the 
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context of magnetogenesis (see first paper of Ref. ra and references therein). A delayed 
reheating has the effect of increasing N max . The largest value of N max in the case of a stiff 
pos-inflationary phase can be estimated as N max = 78.3 + (1/3) lne. In the sudden reheating 
approximation we have iV max ~ 63.25 + 0.25lne which is numerically close to the minimal 
number of efolds N min needed to solve the kinematic problems of the standard cosmological 
model (i.e. N min ~ A max ). Recall, as usual, that An denotes the fiducial amplitude of the 
scalar power spectrum at the pivot scale k p = 0.002 Mpc -1 . We shall not deviate from these 
fiducial values and assume the standard concordance lore for all the late time cosmological 
parameters. 

In Fig. [3]we present the contours of constant log Q E (the common logarithm critical frac¬ 
tion of electric energy density) and the contours of constant logfls (the common logarithm 
of the critical fraction of the magnetic energy density). In Fig. [3] the spectra have been 
evaluated at the maximal wavenumber kr — 1. This is strictly speaking not the highest 
frequency of the electric power spectrum which is exponentially suppressed as soon as the 
reheating start^J To extend artificially the electric power spectrum beyond the frequency 
fixed by the postinflationary conductivity (for more details see [221, 25]) implies that the 
resulting constraints are much stronger. Since we are here trying to constrain the scenario it 
makes sense to impose here the most demanding requirements at the most demanding scale. 


In Fig. [3] the values of the rates encompass all the four quadrants of Figs. [T] and [2] This 
choice has been made also for illustrative reasons since, as explained above, in the remaining 
figures of the section we shall limit our attention to the first quadrant. The results of Fig. 
[3] illustrate the domains of the parameter space where the corresponding critical fractions 
are still perturbative around the maximal wavenumber of the spectra. The contours where 
one of the two critical fractions is larger than 1 should be ideally excised from the physical 
region of the parameter space. We note that the classes of models with ff — 1 seem to 
be comparatively less constrained than the ones with f t = 1. In Fig. [3] F E = F E + 1 and 
F E = 1/2 have been excisec 15 In Fig. [ 3 ] and in the forthcoming figures we have chosen, as 
already mentioned, An = 2.41 x 10~ 9 and e = 0.01. 

The results of Fig. [3] also fix the allowed excursion of F e and F p . Indeed at the maximal 
frequency of the spectrum the amplitude depends also on the maximal value of F E and F p . 
In Fig. [3j for illustrative reasons, we used values \F B \ < 2 and \F E \ < 2. Each value of /) 
(or ff) is compatible with a maximal range of F E and F B . For instance, if ff = 0(1 ) we 
have that, at most \F E \ < 5 and \F E \ < 5: larger values imply that Q E ~ 0(1) around the 
maximal wavenumber of the spectrum. 

14 Imposing the correct boundary conditions between the inflationary phase and the postinflationary phase 
entails also an overall suppression of the electric power spectrum in comparison with the magnetic power 
spectrum [25] . _ _ 

15 As discussed in Eqs. (3.29) and (3.30) this is due to a singularity of the parametrization (not of the 
model). 
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Figure 3: The contours for constant logO B and for constant logO#. The plots on the left 
refer to the case // = 1 and for N t = N mSLX = 62.09; the plots on the right refer to the case 
fi = 1. The various labels illustrate the values of the critical fractions on the given curve. 
On the horizontal axis we illustrate the rate of variation of while on the vertical axis the 
rate of variation of qe is illustrated. 
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5.2 Backreaction constraints 


If the total energy density of the large-scale modes does not exceed the critical energy density 
the following integral 

r k f T dx 


*^(fii d ) 


' kiT X 

2vr 


= —eAn /, 


Q B (a^)x 5 ~ 2 ^ + Q E (a^)x 5 ~ 2 I- 1 ' 

2//(|(j|-l) 


W l-i 


(5.3) 


cannot exceed 1. In practice it will also have to be much smaller than 1, i.e. at most 10~ 3 . 
In the most constraining situation the integral goes from hi = l/r* (corresponding to the 
mode leaving the Hubble radius when a = a*) up to kf — 1/t/ (corresponding to the mode 
reentering the Hubble radius when at a = a,f). Performing explicitly the integrals in the case 

7^ 5, Eq. (5.3[) implies 


2|cr| 7^ 5 and 2|cr — 1 
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5 ~ 2|CT ” 11 ^ a A 5_2|CT_11 
Aa// 


• (5-4) 


To analyze the requirements imposed by Eq. (5.4) we shall first consider the limit Qe —> 1 
and <2# —» 1. With this approximation it will be easier to pin down the area of the allowed 
region in the (F Bl Fe ) plane. This preliminary step will then be supplemented by the 
accurate contour plots obtained with the exact form of Q E and Q B . Indeed the corrections 
due to Qe and Q B can be neglected in the first approximation but they are not irrelevant. 

Thus, in the limit Qe —» 1 and Q B —>■ 1, Eq. (5.4) implies that the following three 
inequalities must be separately satisfied in the physical region of the parameter space: 

< 0 , 

— 5 + 2|cr| < 0, 
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r In fi 1 

U H 


5 + 2|cr — 1| < 0, 


(5.5) 


where N t denotes, as usual, the total number of efolds. Recalling that In (////,) = 2 (F B — 
F E )N t , the two complementary physical situations correspond to fj = 1 and to = 1. 
Indeed these two cases well represent the effect of the initial conditions on the allowed region 
of the parameter space. If fi ^ 1 (but f % ~ 0(1)) the situation is qualitatively close to the 
case fi — 1; if /* 1 the qualitative features of the spectra are close to the case ff — 1. 

16 The cases 2|<r| = 5 and 2|er — 1| = 5 must be separately treated since the integral lead to logarithmic 
divergences. More specifically, if cr = 5/2 and ff = 1, then everything is okay besides the logarithmic growth 
which is not essential for the bound. Conversely, if cr = 5/2 and fi = l we must anyway demand /r < 0. If 
a = —5/2 and ff = 1 the spectra are excluded since they get easily overcritical. In the case a = —5/2 and 
fi = 1 the spectra are viable provided 3/r + 2 < 0. When cr = 7/2 and // = 1 the spectra are excluded while 
for fi = 1 they are only compatible provided + 2 < 0. Finally, if cr = —3/2 the constraints are under 
control for // = 1 (and for f, = 1 but provided fi < 0). 
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Figure 4: The two contour plots illustrate the common logarithm of Ft E without any sim¬ 
plifying approximation. In the left plot the region of Eq. (5.8) is illustrated. The plot on 
the right refers instead to the region of Eq. (5.9). In both cases the wavenumber has been 
chosen to be 10 _4 Mpc _1 . 


5.2.1 The case // = 1 

In this case he backreaction constraints can be evaded provided: 


5 — 2|cr| > 0, 


5-21(7-11 > 0. 


(5.6) 


When o > 1 the conditions imposed by Eq. (5.6) imply, respectively, a < 5/2 and a < 7/2 


Thus, in the domain a > 1 the most constraining bound is a < 5/2. In the region 0 < a < 1 
Eq. (5.6) demands that —3/2 < cr < 5/2. But both conditions are always verified in the 


interval 0 < cr < 1. Therefore no supplementary bounds arise in this case. Finally, if a < 0 


Eq. (5.6) demands cr > —5/2 and a > —3/2. The most constraining of the two previous 
conditions is cr > —3/2. In Fig. [2] the three domains cr>l, 0<cr<l and a < 0 have been 
specifically illustrated in the (F B , F E ) plane. 

Collecting together all the relevant requirements, in the case // = 1 the bounds of Eqs. 


(5.3) and (5.4) imply that —3/2 < cr < 5/2. Recalling Eq. (4.2) we then have the following 
requirement in the (F B , F E ) plane: 


-3 < 


1 - 2F f 


1 + F B - Fj 


< 5. 


(5.7) 


The inequalities of Eq. (5.7) are veriEed in the following two non-overlapping regions of the 
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first quadrant (i.e. F E > 0 and Fb > 0) of the (Fb, Fg) plane: 


Fe > Fb + 1, Fb > -Fb + -, 


Fe < 1 + Fb, 


^ 4^4 


(5.8) 

(5.9) 


Notice that in the region of Eq. (5.8) the electric coupling grows always faster than the 


magnetic coupling (i.e. F E > F B ) and consequently / decreases (i.e. F < 0). Conversely 


in the region determined by Eq. (5.9) we have that F E > F b when 0 < Fb < 2 and 


Fe < F e whenever Fb > 2. The straight line Fe = 5Fb/3 + 4/3 corresponds to the case of 
a flat magnetic power spectrum (i.e. ns —>■ 1). This is an isospectral line in the sense that 
along this line the magnetic and the electric spectral indices do not change. Conversely the 
amplitude of the spectrum still has some dependence on Fb- 


log flfi, ff— 1, N t — N max log fig, ff— 1, N t — N max 



Figure 5: The contours of constant logf^B in the case of Eq. (5.8) (plot at the right) and 


in the case of Eq. (5.9) (plot at the left); the typical scale has been chosen, as in Fig. [fj to 
coincide with k = 1 Mpc^ 1 . 


The results of Eqs. (5.8) and (5.9) are conhrmed by the accurate analysis of the relevant 


contour plots. In the left and right plots of Fig. [4] we illustrate the common logarithm of fl E 


in the case of the regions discussed, respectively, in Eqs. (5.8) and (5.9). The power spectra 


are evaluated for wavenumbers comparable with the infrared branch of the spectrum, i.e. 
typically of the order of 10 -4 Mpc -1 . Larger wavenumbers relevant for the magnetogenesis 
problem (i.e. k = 0(Mpc -1 )) lead to slightly larger figures which are however always of the 
same order. This is due to the quasi-flatness of the spectra in this domain. At the highest 










































frequency the plots of Fig. [3] apply and, even there, the energy densities are much smaller 
than the critical one. In both plots N t = N max = 62.09; the various labels illustrate the 
values of the critical fraction on the given curve. The shaded areas illustrates the regions 


where, according to Eqs. (5.8) and (5.9) the backreaction constraints are satisfied. 


In Fig. [5] we illustrate the same cases discussed in Fig. [4] but in terms of the common 
logarithm of f 1b- The fiducial scale is always k = 10 -4 Mpc -1 and the other parameters 
are also unchanged in comparison with Fig. |4j The combination of Figs. IB and [5] shows 
that the backreaction constraints are satisfied both for the electric and for the magnetic 
fields in the physical range of wavenumbers. In the last part of this section the safe regions 
pinned down by the backreaction constraints will be confronted with the magnetogenesis 
requirements. 


5.2.2 The case f t — 1 


We are now going to address the complementary case where /* = 1. According to Eq. (5.5) 
in the case /« = 1 the following inequalities must be satisfied: 


—2/i(|cr| - 1) + 5 - 2|cr| > 0, 

—2//( rr — 1) — 5 — 2 fx — 1 >0, 
—2/i(|cr| - 1) > 0. 


(5.10) 


If cr > 1, Eq. ( |5.10[ ) implies that the two most constraining inequalities are 2/i(cr — 1) < 0 
and 5 — 2cr > 2/j>{a — 1). In the plane (F B , Fe) this requirement translates into the following 
triplet of inequalities: 


Fe > F b + 1 , Fe < Fb , Fb > — 2 ; 


(5.11) 


notice that F b > —2 is automatically satisfied since we work in the region F b > 0 and 

Fe > 0. 


If a < 0 the conditions of Eq. (15.10 ) become 


2+ 1) > 0, 5 + 2 a > —2 fi(a + 1), 5 + 2cr > 2 — 2/r(cr + 1). 


(5.12) 


The last inequality of Eq. (5.12) is the most constraining so that —3/2 — /i(cr + 1) < a < 0. 
In terms of Fe and F b this range of cr translates into: 


1 - 2 Fp 


< 0 , 


{Fb ~ Fe ){3 + 2 Fb — 4i+ 


> 0 , 


FEFl+^-f- (5.13) 


2(1 + Fb — Fe) ’ 1 + Fb — Fe 

If 1 + Fb — Fe > 0, then Fe >1/2 and the inequalities defining the allowed regions arc 1 ' 

{Fb — Fe ){3 + 2 Fb — 4 Fe) > 0 , 


Fe < 1 + Fb/ 2. 


(5.14) 


17 If 1 — F Fj + F b < 0, Eq. (5.13) demands that F E < 1/2 which is not compatible with F E > F B + 1. So 
this possibility must be discarded. 
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Figure 6: The common logarithm of tt E hi the region (5.15) (plot on the left) and in the 
region (5.16) (plot on the right). The same parameters and the same scales of Figs. [4] and 
[5] have been employed. 

For Fb > Fe the allowed region is: 

1/2 < F b < 3/2, 1/2 < F E < F B , 


For F b < F e the allowed region is: 


Finally in the region 0 < a < 1 the backreaction constraints are satisfied in the region 
F e < 1/2. 

In Figs. [6]and[7]we illustrate, respectively, the common logarithm of Vt E and of Vt E - The 


3/2 <F e < 3/4 + F b /2. 

(5.15) 

3/4 + F b /2 < F e < 1 + F B /2 

(5.16) 


contour plots on the left refer to the region of Eq. (5.15) while the contour plots on the right 


illustrate the parameter range of Eq. (5.16). The fiducial set of cosmological parameters is 
the same as in the previous figures. 


5.3 Dualizing the contour plots 

The region of the parameter space explored so far coincides with the first quadrant of the 
(F e , F e ) plane. Using duality the contours plots of the other portions of the parameter space 
can be easily obtained and studied, as specifically discussed in sec. [4} Let us just give an 
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Figure 7: The common logarithm of fig in the region (5.15) (plot on the left) and in the 
region (5.16) (plot on the right). The same parameters and the same scales of Figs. [4] and 
[5] have been employed. 


example of this strategy. Let us suppose, for instance, to be interested in the case where 
both rates are negatives, namely Fb < 0 and Fe < 0, as it happens in the third quadrant of 

Figs. [T| and [2j 

Recalling therefore Eqs. (4.11) and ( 4. 12| ) (and also Tabs. [I] and [2]) we can easily dualize 
all the contour plots obtained so far in the regions where the backreaction constraints are 
safely satisfied. In Figs. [8]abd[9]we illustrate the contours plots dual to the ones reported 
in Figs. [4] and [5} Since, under duality, Fb —> —Fe and Fe —> —Fb the axes of Figs. [8] an d@ 
are interchanged in comparison with the axes of Figs. [4] and [5j The same exercise discussed 
in the case // = 1 can be also carried on in all the other relevant cases of phenomenological 
interest and also in the remaining quadrants of the parameter space. 

The plots reported in Figs. [8] and [9] can be related, by symmetry, to the plots of Figs. [4] 
and [5] Of course also the domains of the spectrum must be dualized as specifically discussed 
in section [4j It is however appropriate to stress that Figs. [8] and [9] have been obtained by 
plotting explicitly the relevant critical fractions in the dual regions of Figs. [4] and [5] The 
obtained results show indeed that, a posteriori, duality can be directly applied also to the 
exclusion plots without plotting explicitly the relevant spectra. 
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Figure 8 : We illustrate the contour plots obtained by duality from Fig. [4j Note that the 
axes are interchanged thanks to duality. The notations and the parameters are the same as 

in Fig. [|J 


log H 5 , ff— 1, N t - N max log ff— 1, N t —N max 



Figure 9: We illustrate the contour plots obtained by duality from Fig. [5] Note that the 
axes are interchanged thanks to duality.The notations and the parameters are the same as 
in Fig. [5] 
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5.4 Magnetogenesis requirements 

Having determined the regions of the parameter space where the backreaction constraints 
are enforced we can address the problem of the magnetogenesis constraints. Let us therefore 
focus our attention on the first quadrant, as in the previous part of this section. In units of 
nG 2 the magnetic power spectrum can be written as follows: 


Pb 

nG 2 
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10 


-3.05 


h 2 0 n 


RO 


4.15 x 10~ 5 / V2.41 x 10” 9 


An 


0.01 


Qb^B, h) 





p Nmaxi^B 1) 


y(4—n B )/2 ^ fiNt(4:—n B ) 


(5.17) 


The magnetogenesis requirements roughly demand that the magnetic fields at the time of the 
gravitational collapse of the protogalaxy should be approximately larger than a (minimal) 
power spectrum which can be estimated between 10~ 32 nG and 10 -22 nG. The most optimistic 
estimate is derived by assuming that every rotation of the galaxy would increase the magnetic 
held of one efold. The number of galactic rotations since the collapse of the protogalaxy can 
be estimated between 30 and 35, leading approximately to a purported growth of 13 orders 
of magnitude. During collapse of the protogalaxy compressional amplification will increase 
the held of about 5 orders of magnitude. Thus the required seed held at the onset of the 
gravitational collapse must be, at least, as large as 10~ 15 nG or, more realistically, larger 
than 10~ n nG. 

The simplest way to implement the magnetogenesis requirements is therefore to plot 
the magnetic power spectrum in units of nG 2 for all the regions where the backreaction 
constraints are satished. The areas of the parameter space where 



22 < A < 32, 


(5.18) 


will therefore offer viable models of magnetogenesis. Depending on the sign of the rates and 
on their relative magnitude it will therefore be possible to see whether the selected models 
will evolve either from a strongly coupled or from a weakly coupled regime. 

Let us therefore start with the case where f/ = 1. The relevant regions of the parameter 
space have been already discussed in Figs. [4] and [5j For the same regions we shall now plot 
the common logarithm of Pb/iiG 2 and check if and where the condition of Eq. (5.18) is 
satished. The result of this procedure is reported in Fig. [10] where the labels in the contour 
plot refer to the values of the power spectrum in units of nG 2 . The plot on the left in Fig. 
10 refers to the same region discussed in the plots on the left in Figs. [4] and [5] similarly the 


plots on the right correspond to the plots on the right in Figs. [4] and [5j We clearly see that 
the requirement of Eq. (5.18) is satished in two distinct regions: a stripe of the left plot and 
a smaller stripe in the plot on the right. According to Fig. [4| however, the region in the plot 
on the right must be excluded since in that part the electric helds are overcritical. We are 
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Figure 10: We illustrate the magnetic power spectrum in units of nG 2 for the regions dis¬ 
cussed in Figs. [4] and [5] where the backreaction constraints are all satisfied. 


therefore left with a region which can be analytically expressed as 


5F B /3 + 4/3 < F E < 
)C(X,k,N max J f ) = 


1C 


6 f b + k 


jC-A a ' /C — 4 ’ 
3.05 — A + log ( k/H q) 


Nmax log e - 2 log ff 


log (k/H 0 ) - N max log e - 0.5 log f f 


(5.19) 

(5.20) 


In the case of the fiducial set of parameters corresponding to the previous figures we have 
that the above slice becomes 


5F B /3 + 4/3 <F e < 1.46 + 1.9LF B , A = 22, (5.21) 

5Fb/3 + 4/3 < F e < 1.56 + 2.13Fb, A = 32. (5.22) 


The lower bound in Eq. (5.19) corresponds to the slightly distorted straight line appearing 
both in Figs. [4] and [TO} The upper bound has been obtained from the approximate shape of 
the slice defined by the corresponding contour plot. In this region both gauge couplings are 
increasing (no strong coupling problem) and furthermore / decreases from an initially large 
value to fj = 1. 

The same analysis discussed in the case ff — 1 can now be extended to the case /* = 1 
already discussed in Fig. [6] and [7} The result of this analysis is reported in Fig. 11 where 
we clearly see that the magnetogenesis constraints are not satisfied. 

As already mentioned in the case of the backreaction constraints, the various regions of 
the parameter space can be easily dualized with the purpose of discussing a given region of 
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Figure 11: We illustrate the magnetic power spectrum in units of nG 2 for the regions dis¬ 
cussed in Figs. [6] and [7] where the backreaction constraints are all satisfied. 


fc? 


log /VnG 2 , //= 1, N,= N max 


log P B /nG 2 , f f = 1, N t =N ma]t 



Figure 12: We illustrate the magnetic power spectrum in units of 11 G 2 for the regions dis¬ 
cussed in Figs. [8] and [9] where the backreaction constraints are all satisfied. 

the parameter space once the results on the first quadrant are known. An example along 
this line is obtained in Fig. [12] where the magnetogenesis constraints are illustrated in the 
case of the dual regions already studied in Figs. [8] and [9j We can therefore see that also 
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the region of Fig. 12 leads to a viable class of magnetogenesis models where, however, the 
rates are negative and the gauge couplings decrease rather than increasing. There are some 
who say that these initial data are unnatural. Some other authors claim instead that the 
gauge coupling should follow the gravitational coupling. In conventional inflationary models 
the strong coupling is at the beginning (since the singularity is at the beginning) while in 
bouncing models the singularity is at the end of the inflationary phase. So the most natural 
choice would be to ask for decreasing gauge couplings in conventional inflationary models 
and increasing gauge couplings in bouncing models (see, e.g. third paper in Ref. [25]). 
In this scenario both options can be realized in non-overlapping domains of the parameter 
space. 

The cases ff — 1 and fi — 1 represent in practice the two extremes of the whole range 
of physical models where fi > 1. The cases where f, = 0( 1) are very close to the patterns 
discussed in the /) = 1 case. Conversely when /) » 1 the situation is similar to the case 
ff = 0( 1). It is clear that by appropriately scanning the parameter space we could easily pass 
from an allowed region in the plane (F B , Fe) to a tuning volume in the space (Fb, Fe, fi)- 
In that space the range of variation will be given exactly by the excursion discussed in the 
present paper, i.e. 1 < /,; < where f/ max) coincides with exp [(Fe — F B )N t \ and it 

corresponds to ff = 1. 

Before closing this section let us also comment on the possible excursion of N max . As 
already mentioned if the reheating is delayed N max can increase and get larger by even 15 
efolds [25i , 26]. In this case the effect on the allowed region is minimal. For instance if 
N m ax = 78 the area of Eq. (5.21) becomes Fb /3 + 4/3 < F E < 1.42 + 1.84Fb which is not 
different, for all practical purposes, from the previous expression. 


6 Concluding remarks 

Magnetogenesis models based on derivative couplings have been comprehensively analyzed. 
A similar kind of framework arises in the relativistic generalization of Van der Waals in¬ 
teractions. The key aspect of this class of models is that the electric and the magnetic 
susceptibilities are not bound to coincide all the time, i.e. xe Xb- This implies that also 
the corresponding gauge couplings g E and g B may evolve in time at different rates. After 
presenting a general decomposition of the various coupling functions, the quantization of 
the system has been approached in terms of a newly defined time coordinate reducing to 
conformal time in the case when the electric and the magnetic susceptibilities coincide. 

It turns out that the power spectra depend on three different quantities: the normalized 
rates of variation of the electric and magnetic gauge couplings (i.e. F E and F B ) and also 
fi, i.e the amplitude of g B /g E at the onset of the dynamical evolution that has been taken 
to coincide with the beginning of the inflationary phase. A fourth parameter measures the 
strength of one of the two gauge couplings at the beginning of inflation. The parameter space 
of the model has been scanned by using duality in order to relate the different quadrants 
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of the (F B , Fe ) plane. The physical content of the whole parameter space can therefore be 
reduced, via duality, to the careful analysis of the region where Fe and Fb are both positive 
semidehnite. 

While the arbitrary variation of Fb and Fe is probably the most relevant aspect, a full 
account of the arbitrary variation of f t is not particularly significant at this stage of the 
analysis of the model. To explain the general trend it is sufficient to examine the case where 
fi = 0( 1) and the case where // = 0(1), recalling that // denotes the value of / at the end 
of inflation. The parameter space of this magnetogenesis scenario has then be accurately 
charted in the (Fb, Fe) plane. As a consequence of this analysis the models where fi = 0(1) 
are favoured when at least one of the gauge couplings is initially strong while the models 
ff = 0(1) are favoured when both gauge couplings are small at the beginning of inflation. 

There exist wide regions in the parameter space where backreaction effects are negligible, 
gauge couplings are initially small and ff = 0(1). This area is defined by the corresponding 
excursion of the normalized rate^jin the (Fb, Fe) plane: 

5F b /3 + 4/3 < F e < 1.46 + 1.91+b, 0<F B <5. (6.1) 


The functional dependence of the magnetic spectral index upon F b and Fe changes from 


quadrant to quadrant (and also within the same quadrant). In the area defined by Eq. (6.1) 
it is given by: 

5 + 6 F b - 4:F e f ^ 

kb = -i,p _ t? » 1 <n B < 1.9. ( 6 . 2 ) 

-L i - 1 fj 1 fff 

The second relation appearing in Eq. (6.2) follows from Eq. (6.1) and from the explicit 


expression of n B in the allowed region. In the range of Eqs. (6.1) and (6.2) the magnetic 
power spectrum at the scale of the protogalactic collapse is always larger than 10~ 22 nG 2 . 
As in the case of curvature perturbations, the scale-invariant limit of the magnetic power 
spectrum holds for n B —> 1 where its amplitude is: 


P ; 


B 


lim 

1 nG‘ 


= 10 


-2.805 


h 2 0 n 


R0 


An 


4.15 x 10-V V2.41 x 10” 9 / V0.01 


1 + 2 F, 


B 


(6.3) 


The scale-invariant amplitude still depends on F b . If the magnetogenesis constraints are 
relaxed by requiring that the magnetic power spectrum at the protogalactic collapse exceeds 
only 10 -32 nG 2 the range of allowed values of n B gets slightly wider, namely 1 < n B < 2.2. 

In Tab. [3] we illustrate the values of the magnetic power spectra in the allowed region 
defined by Eqs. (6.1) and (6.2). In Tab. [dj for the same portion of the parameter space, we 
report the magnetic spectral indices. Notice that the boundary lines of the allowed region 
are isospectral in the sense that, along these lines, the magnetic spectral index does not 
change or it changes very little. It is interesting to stress that, in this model, ff 1 implies 
the equality of the gauge couplings at the end of inflation but not necessarily the coincidence 


ls The curves defining the allowed region depend on the total number of efolds which will be taken to 
coincide with iV max ~ 63.25 + 0.25 In e , i.e. the maximal number of efolds today accessible by observations 


37 

















F b 

F e = 5F b /3 + 4/3 

F e = 1.46 + 1.91Fb 

F e = 1.56 + 2.13+b 

0 

s/Pn = lO" 2 - 80 nG 

yfPs = lO" 12 - 25 nG 

y/P^ = io~ 16 - 69 nG 

2 

y/FE = 10- L40 nG 

\[Pb = lO" 1091 nG 

\[Pb = lO -15,86 nG 

4 

\/Pb — 10 -0 ' 89 nG 

s/Pb = lO" 10 - 48 nG 

JJT = 10 -15.52 llG 


Table 3: Typical values of the magnetic power spectra in the allowed region and at the 
fiducial scale k = lMpc -1 . 


F b 

Fe — 5Fb/3 + 4/3 

F e = 1.46 + 1.91Fb 

F e = 1.56 + 2.13+b 

0 

n B = 1 

tib = 1.826 

tib = 2.214 

2 

n B = 1 

77 B = 1.807 

n B = 2.226 

4 

riB = 1 

77 B — 1.804 

tib = 2.228 


Table 4: Typical values of the magnetic spectral index in the allowed region. 


of the rates. On the contrary, if the rates are coincident during inflation (i.e. Fe —y F b ) 
the conventional situation is recovered and the two gauge couplings evolve at the same rate. 
The curve Fe = Fb is not isospectral and it is not included in the area defined by Eqs. (6.1) 
and (6.2). This occurrence shows, once more, that when F E ^ F B the parameter space gets 
indeed wider in comparison with the conventional situation. 

All in all if the magnetic and the electric susceptibilities do not coincide the allowed 
regions in the parameter space of inflationary magnetogenesis gets wider in comparison 
with the conventional class of models. According to some, the natural initial conditions 
for inflationary magnetogenesis strictly demand minute gauge couplings at the beginning of 
inflation and larger gauge couplings at reheating. According to a different way of thinking 
gauge coupling should follow the gravitational coupling: as the production of a flat spectrum 
of curvature perturbations demands a strong gravitational coupling in the past, similarly a 
quasiflat magnetic field spectrum is realized in the case of a decreasing gauge coupling which 
gets progressively smaller during inflation. In this scenario both options are plausible in 
different regions of the parameter space. 
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